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The heights of iterates of the discrete Painleve equations over number fields appear to grow no 
faster than polynomials while the heights of generic solutions of non-integrable discrete equations 
grow exponentially. This gives rise to a simple and effective numerical test for the integrability 
of discrete equations. Numerical evidence and theoretical results are presented. Connections with 
other tests for integrability and Vojta's dictionary are discussed. 
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Over the past decade and a half several criteria have been suggested as detectors of integrability for maps and 
discrete equations. Many of these criteria echo the observation of Veselov that . . . integrability has an essential 
correlation with the weak growth of certain characteristics Q]. A number of authors have studied rational maps and 

\ discrete equations for which the degree of the n^ 1 iterate y n as a rational function of the initial conditions grows no 
faster than a polynomial in n (JQ-Q In particular, the algebraic entropy introduced by Hietarinta and Viallet Q is 
a measure of this degree growth and is related to Arnold's idea of complexity. 

The singularity confinement property of Grammaticos, Ramani, and Papageorgiou has led to the discovery of many 
i— i" integrable discrete equations Q. Hietarinta and Viallet have shown that there are non-integrable equations that 
possess the singularity confinement property, which led them to add the (growth-type) condition of zero algebraic 
. • entropy. 

It was suggested in [j| that the existence of sufficiently many finite-order meromorphic solutions of a difference 
r— h ■ equation is a natural analogue of the Painleve property and a detector of integrability. It has been shown in 6] that 

1 " , '_ if an equation of the form y(z + 1) +y(z — 1) = R(z, y(z)), where R is rational in its arguments, admits a non-rational 

finite-order meromorphic solution then either the equation can be transformed to one of the known discrete Painleve 
t-H ' equations or y also satisfies a (first-order) discrete Riccati equation. 

Important connections between the differential and discrete Painleve equations, representations of affine Weyl 
groups, and the geometry of certain rational surfaces have been found in 7]. Costin and Kruskal have suggested that 
the theory of analyzable functions is the appropriate language in which to describe the Painleve property for discrete 
equations @ . Roberts and Vivaldi have considered maps over finite fields and used orbit statistics to single out detect 
\ maps with a polynomial integral of motion [j| . 

■ In the present letter a slow-growth property is described which is very easy to test numerically. It involves considering 
the iterates of a discrete equation in an appropriate number field (i.e., a finite extension of the rationals) and examining 
the growth of the height of these iterates. The height H(x) of an element a; of a number field k is a measure of the 
• complexity of x. 

We will, for the most part, only deal with the case k = Q. The height of a non-zero rational number x € Q is 
H(x) — max{|p|, |g|}, where x — p/q and p and q have no common factors. The height of is defined to be H(0) = 1. 
We begin by considering the growth of heights of iterates y n of the equation 



y n+1 + y„_i = — + b n , (1) 

d 1 Vn 
_ i 

where j/o and yi are given rational numbers and a n and b n are chosen to be in Q for all n£Z. This guarantees that all 
finite iterates y n are also rational numbers. In FIG.^ l°gl°g-ff(yn) has been plotted against logn for three solutions 
of equation JJJ. In each case, the initial conditions are yo = 2/5, y\ = 3/7 but the choices of a n and b n differ. 

When a n and b n are constants, equation can be solved in terms of elliptic functions. When a n = An + ji and 
b n = v, for constants A, fi and v, then equation Q is an integrable discrete equation related to the first Painleve 
equation. In FIG.^ there are two integrable cases (a„ = 3, b n = 5 and a n = 3 + n, b n = 5) corresponding to 
the asymptotically straight line plots, while the third non-integrable case (a n = 3, b n = 5 + n) corresponds to the 
asymptotically non-linear curve. 

The above example motivates defining a polynomial discrete equation such as Q to be Diophantine integrable if 
the logarithmic height of iterates, h(y n ) = log H(y n ), grows no faster than a polynomial in n. This idea is certainly 
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FIG. 1: Plot of loglogif(j/ n ) vs logn for equation Q 



related (and possibly equivalent) to the degree growth/algebraic entropy approaches described above but is much 
quicker to check numerically for a large number of iterates. For rational maps over k = Q, Abarenkova, Angles 
d'Auriac, Boukraa, Hassani and Maillard 01 used the heights of iterates as a measure of complexity and to find 
special values of a parameter for which a particular map was integrable. Below we will see that for general rational 
discrete equations we are led to the growth of heights over arbitrary number fields from an analogy with Nevanlinna 
theory. We will use this analogy to prove that a large class of discrete equations are not Diophantine integrable. 
However, in contrast to the existence of finite-order meromorphic solutions as a detector of integrability for difference 
equations, it is very simple and quick to use a symbolic computing package to test (numerically) for Diophantine 
integrability. 

To deal with the case in which an iterate becomes infinite it is natural to work in projective space, however, for 
the purposes of this letter we will only consider finite iterates. Many autonomous versions of the discrete Painleve 
equations are known to be solved in terms of elliptic functions. In fact, these equations are essentially the addition 
law on the cubic. For any infinite sequence of rational points on an elliptic curve such that any iterate is the sum (on 
the cubic) of the previous two, the logarithmic height grows like n 2 . 

Next we consider the so-called qPvi equation, which is the system 

fnfn+l = 9n+l ~ Ctq n+1 gn+1 - /?g" +1 ,^ 

cd g n+ i - 7 g n+ i - S 

9n9n+l fn - aq n f n - bq n 



jS f n - c f n - d 

subject to the constraint q = af3~/5 / abed. The system was discovered by Jimbo and Sakai as the compatibility 
condition for an isomonodromy problem llj and is an integrable discretization of the sixth Painleve equation (Pvi)- 
FIG. is a plot of log log max{ #(/„), H(g n )} against logn for iterates of equation © with the initial conditions 
fo = 2/3, g = 3/4 and the choice of parameters (a, /3, 7, S, a, b, c, d) = (15/7,4/3,1/2,1,8/7,5/7,2,1/7). The two 
graphs represent two different choices for q, namely, q — l/2(= af3^S / abed, i.e., the integrable case corresponding to 
the asymptotically linear graph) and q — 2. Once again we see that the logarithmic heights h(f n ) and h(g n ) appear 
to grow polynomially in the integrable case and exponentially in the non-integrable case. 

Now we will discuss some fundamental identities concerning heights and discrete equations. Let 

a + aix + ■ ■ • + a p x p 



6 + bix + • • • + b q x q ' 
be an irreducible rational function of x of degree d — max{p, q}. Then 

dH{x) d < H(R) < C 2 H(x) d , (3) 



FIG. 2: Plot of log log max{H(f n ), H(g n )} vs logn for equation 



where C\ and C2 are polynomials in the heights of the coefficients a{, bj. The second inequality in @ is straightforward. 
In the first inequality C\ is proportional to the resultant of the denominator and numerator of R. So the logarithmic 
height h{ ■ ) = \ogH{ ■ ) satisfies 

\h{R) -dh(x)\ <logC, (4) 

where C is a polynomial in H{a,i) and H(bj). 
Consider the first-order discrete equation 

_ ao(n) + ai(n)y n + h a P (n)yP 

n+ b (n) + bi(n)y n H h&g(n)j/&' 

where the a^'s and bj's are in Q[n] and a p and 6 g are both not identically zero. It follows that for all integers n larger 
than some no, the degree of the right side of equation JSJ as a function of y n is independent of n. If the right side 
of equation J5j is written in irreducible form for n > no then this degree is d — max{p, q}. Taking the logarithmic 
height of both sides of equation JSJ and using equation Q gives 

h{Vn+i) = dh(y„) + O(logn). 

So if H(y n ) grows faster than any polynomial in n (i.e., iih(y n ) grows faster than logn) then h(y n ) grows exponentially 
unless d < 1. In this case we are left with the (integrable) discrete Riccati equation y n +i = {ao(n)+ai(n)y n }/{bo(n) + 
bi(n)y n }, which can be solved via a second-order linear discrete equation. Hence the demand that solutions grow 
no faster than polynomials singles out the discrete Riccati equation in the same way that the Painleve property 
for differential equations singles out the differential Riccati equation. Note that any periodic orbit of equation JSJ), 
% n ) = 0(lV 

Osgood ,12$ observed that there is an uncanny formal similarity between the basic definitions and theorems of 
Diophantine approximation and those of Nevanlinna theory. Independently, Vojta [T^ | constructed a dictionary which 
provides a detailed heuristic for the "translation" of concepts and propositions between the two theories. 

The fundamental idea of Nevanlinna theory is that much information about a meromorphic function / is obtained 
by studying a kind of averaged behaviour of / on the disc D r := {z : \z\ < r}. The main tool is the Nevanlinna 
characteristic T(r, /). For an entire function, T{r, /) behaves like logmax| z | =r \f(z)\. In general, T(r, /) is the sum 
of two terms, one of which is the average of a certain function of |/| on \z\ ~ r and the second is a measure of the 
number of poles of / in the disc D r . 

According to Vojta's dictionary, a statement in Nevanlinna theory about the Nevanlinna characteristic of a mero- 
morphic function corresponds to a statement in Diophantine approximation about an infinite set of numbers in a 
number field. In |5( it has been argued that a natural analogue of the Painleve property for difference equations is the 
existence of sufficiently many finite-order meromorphic solutions. Suppose that, in going from a difference equation 
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for y{z) (in the complex plane) to the corresponding discrete equation for y n (in which the independent variable 
is restricted to the integers), we find that all iterates y n are in some number field k for initial values chosen in k. 
Via Vojta's dictionary, the statement that y(z) is a finite-order meromorphic function corresponds to the statement 
that h(y n ) grows no faster than a polynomial. Hence Diophantinc intcgrability is the natural analogue for discrete 
equations of the finite-order growth Painleve- type condition for difference equations. 

We now derive a number of results about Diophantine integrability by exploiting this formal similarity with 
Nevanlinna theory. It is straightforward to check that for any rational numbers xi, . . . ,xn, the height satisfies 

H (^2f=i x j) — NYljLi H(xj) and H M^Li x j*j ^ Y[f=i-^( x j)- ^ follows that the logarithmic height h has the 
following properties. 

(N \ N 
J2 X A ^ X>te)+logW, (6) 

(N \ N 
3=1 / 3=1 

Note that if we replace the rational numbers Xj and the logarithmic height function h( ■ ) in © and JJJ by meromorphic 
functions fj and the Nevanlinna characteristic T(r, ■ ) respectively, we obtain two standard identities in Nevanlinna 
theory. Similarly, it is natural to consider equation to be the Dio pha ntine analogue of the Valiron-Mohonko 
Theorem [l4| . These are some of the fundamental results used in 0, ITU Il6t Il7| to find necessary conditions for the 
existence of finite-order meromorphic solutions of difference equations. 
Consider either of the two second-order discrete equations 

Vn+i +y n -x = R{n,y„) or y n +iy n -i = R(n,y n ), (8) 

where R(n, y n ) is as in equation (JjJJ. Taking the logarithmic height of either of the equations in JH|) and using equation 
(@J and the inequality © or Q yields 

h(y n+ i) + h(y n -i) > dh(y n ) + O(logra). 

It follows that, provided h(y n ) grows faster than O(logn), either it grows faster than any polynomial in n or d := 
deg Vn (R(n, y n )) < 2. The condition d < 2 is consistent with a number of known integrable discrete Painleve equations 
such as the special cases of equation described above. It is also consistent with the equation 



Q» + b n y n 

Vn+l + Vn-l = — j 5—, (9) 

1 ~Vn 



which is integrable when a n is a constant and b„ is linear in n, when it is known as the so-called discrete Painleve II 
equation, dPn. 

FIG. |21 consists of two sequences corresponding to solutions of equation @ in the non- integrable case a n = 2n — 1, 
b n — 2n 2 — 2n + 3. The initial conditions (yo, yx) are (11/12, 11/23) and (11/23, 11/12). The second initial condition 
corresponds to the asymptotically linear case in FIG.|3| This solution is not generic in that it is also a solution of the 
first-order (Riccati) discrete equation y n +i = (n 2 + y n )/( 1 — y n ). 

The conclusion that d < 2 in equations JBJ if y n has "slow growth" (more precisely, h(y n ) grows faster than O(logn) 
but is bounded above by a power of n) is analogous, both in the conclusion and in the proof, to a result proved in Q. 
A number of generalisations of the result in [j| have appeared in 0, . Many of these results also have Diophantine 
analogues, including the results in [l^ related to qPyi (equation [21 above.) 

In order to find Diophantine analogues of all the proofs in ^(|, we need to have one extra identity. Fix N, let I be 
the set I = (1, . . . , N) and for each non-empty JC / assume aj S Q. Then 



It follows that 



h ( e «■> n x i p i j i n H ^ a ^ ) 1 n h ^ 

j ci jeJ ) \jci 



h ( E a J II x i ) ^ E h ^ + 10 S( C )' ( 10 ) 
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FIG. 3: Plot of log log H(y n ) vs logn for equation @ 
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FIG. 4: Plot of log log H{y n ,n) vs logn for equation (II IB 



where C = J2jci h( a j) + l°g l-^li which is independent of the x^s. Inequality (|10|l is a natural Diophantine analogue 
of a useful inequality in Grammaticos, Tamizhmani, Ramani, and Tamizhmani |16| (number 2.11), which is used 
to obtain the general form of many integrable discrete equations. The details of this and other calculations will be 
published elsewhere. 

Finally we consider the lattice equation 



1 



£/m,n+l 2/m+l,n 



(ii) 



where a is a constant. In FIG.0] log log if (y n ,„) has been plotted against logn for two different values of a (a = 1 
is the asymptotically linear case, the other is a — 2) corresponding to the same initial conditions. The case a = 1 is 
known to be integrable while a = 2 is not. 
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